Stretching weakly bending filaments with spontaneous curvature in two dimensions 
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Some important biomolecules (for instance, bacterial FtsZ and eukaryotic DNA) are known to 
posses spontaneous (intrinsic) curvature. Using a simple extension of the wormlike chain model, we 
study the response of a weakly bending filament in two dimensions to a pulling force applied at its 
ends (a configuration common in classical in-vitro experiments and relevant to several in-vivo cell 
cases). The spontaneous curvature of such a chain or filament can in general be arc-length dependent 
and we study a case of sinusoidal variation, from which an arbitrary case can be reconstructed via 
Fourier transformation. We obtain analytic results for the force-extension relationship and the width 
of transverse fluctuations. We show that spontaneous-curvature undulations can affect the force- 
extension behavior even in relatively flexible filaments with a persistence length smaller than the 
contour length. 

PACS numbers: 36.20.Ey,87.15.ad,82.37.Rs 



I. INTRODUCTION 

Semiflexible polymers are fluctuating macromolecules 
with a bending stiffness intermediate between that of 
a random coil (Gaussian chain) and a rigid rod. They 
have been the subject of intense experimental and the- 
oretical study in recent years because many important 
biomolecules fall in this category. The elastic behavior 
of the structural elements of the cytoskeleton (F-actin, 
microtubules, intermediate filaments) as well as that of 
DNA strongly depend on their bending stiffness [TJ [2] . 

A widely used minimal model of semiflexible polymers 
is the wormlike chain model (WLC) which treats the fil- 
ament as a locally inextensible one-dimensional fluctuat- 
ing curve with bending stiffness O 0]. This model in- 
volves only two length scales, namely, the contour length 
L and the persistence length L p , which is the correla- 
tion length of the polymer directedness related to the 
chain bending modulus n = ^L p ksT(d — 1), where d is 
the dimensionality of the embedding space. In the limit 
L/L p ^$> 1 the semiflexible chain transforms into a Gaus- 
sian coil and the question of excluded- volume interaction 
may become relevant; in the opposite limit L/L p <C 1 
the chain approaches the limit of rigid rod, in which 
case the additional modulus penalizing the local stretch- 
ing/compression may come into play [5j. Despite its sim- 
plicity, the WLC model has been successfully used to de- 
scribe some important experimental results such as the 
strong stretching of ds-DNA [6 . However, the advent of 
single molecule experiments has revealed many occasions 
where the complexity of the microscopic architecture of 
the filament appears to affect its macroscopic behavior 
in a way that cannot be accounted for by the classical 
WLC model. 

An important property of some biomolecules which 
goes beyond the basic WLC model is spontaneous lo- 
cal curvature. FtsZ is a protein of the prokaryotic cy- 
toskeleton which plays an essential role in cell division 



as it assembles into the contractile ring (Z-ring) that di- 
vides the bacterium j7j. It can be visualized on a mica 
surface using atomic-force microscopy (AFM). An impor- 
tant property of many FtsZ protofilaments is the spon- 
taneous curvature of their relaxed state. Their role in 
force generation strongly depends on the determination 
of their persistence length [8^ . Assuming the relaxed state 
of FtsZ protofilaments to be straight, Dajkovic et al. [5] 
estimated L p ss 200 nm which is much too small. On the 
other hand, Horger et al. [TU] taking into account the 
finite spontaneous curvature measured by AFM (average 
radius of curvature w 100 nm), determined a persistence 
length L p w 4 fim which is close to the value calculated 
from models of the flexural rigidity. 

The existence of stable bent configurations of DNA has 
been known for many years and they play important roles 
in many basic genetic processes [llj . Spontaneous curva- 
ture in the local structure of DNA is sequence-dependent 
|12) . As in the case of FtsZ, the presence (or absence) of 
local spontaneous curvature in a DNA filament is crucial 
in determining the persistence length and its relation to 
the bending rigidity. In a recent experiment, AFM imag- 
ing of two-dimensional human DNA chains and compar- 
ison to measurements on both intrinsically straight and 
uncorrelated HCV DNA has demonstrated their sponta- 
neous curvature |13j . 

The theoretical study of small (Gaussian) thermal fluc- 
tuations of three-dimensional filaments with spontaneous 
curvature and torsion has yielded analytical expressions 
for the orientational correlation functions and the per- 
sistence length [14] . An auxiliary field theory of three- 
dimensional polymers with intrinsic curvature has ob- 
tained the orientational correlations, mean square sep- 
aration and the force-extension relation of such chains in 
the L/Lp ^ 1 limit, showing that helical configurations 
arise due to the excluded- volume interaction [15]. Real 
filaments with spontaneous curvature in three dimensions 
also possess a finite twist rigidity which has been shown 
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to affect the end-to-end distribution [TOj and the response 
to a stretching force giving rise to instabilities [17l [18] . 
The zero-temperature conformations of two-dimensional 
filaments with constant spontaneous curvature under an 
applied force for various boundary conditions, as well as 
the linear response at finite temperature to an applied 
force, have been recently investigated in [T§] . 

In this paper, we study the response of a two- 
dimensional WLC with spontaneous curvature to a 
stretching force applied at the ends. Although in gen- 
eral one of course expects the stretched chain to explore 
fluctuations in both perpendicular directions, the two- 
dimensional limit is a meaningful simplifying approxima- 
tion that permits us to carry out a complete analytical 
study, while retaining some of the relevant qualitative 
features of the general problem. The three-dimensional 
problem is more complicated in two nontrivial aspects. 
Curvature in 3d is always positive definite whereas in 
2d we can define the signed curvature. In addition, the 
presence of twist rigidity complicates matters even fur- 
ther [TMTB] . As we explain in Sec. VII, however, the 
main conclusion of this work which is the softening of 
the stretching response due to a spatially varying spon- 
taneous curvature is expected to carry through in 3d. Be- 
sides, many imaging experimental investigations |10[ 113] 
are carried out on molecules confined to a surface. A 
stretching experiment on DNA electrostatically confined 
to the surface of a fluid membrane is presented in |20j . 
The elasticity of semiflexible polymers without sponta- 
neous curvature in two dimensions has been studied an- 
alytically in |21j . The spontaneous curvature can in gen- 
eral be inhomogeneous, dependent on the location along 
the arc length s along the chain. Any such variation can 
be decomposed into Fourier harmonics and so we study 
a case of local curvature varying in an arbitrary sinu- 
soidal fashion along the chain length. We assume that 
the filament is weakly bending. This approximation can 
be realized in the case of filaments with small sponta- 
neous curvature and small thermal fluctuations or in the 
case of strong stretching forces which have flattened out 
the initially large spontaneous curvature and/or the large 
thermal undulations. This is the simplest possible exten- 
sion of the WLC model to include spontaneous curvature 
analytically. A formally similar approach has recently 
been employed in [25] to study uncorrelated quenched 
disorder in the spontaneous curvature of stretched two- 
dimensional filaments. 

The paper is organized as follows. In Sec. II we in- 
troduce our model and derive the Gaussian energy func- 
tional which characterizes the stretched polymer. In the 
absence of thermal fluctuations, the polymer assumes the 
ground-state conformation which minimizes the elastic 
energy for given boundary conditions; finding this confor- 
mation is the subject of Sec. III. In Sec. IV we show that 
thermal fluctuations decouple from the athermal mechan- 
ical response and are independent of the spontaneous 




FIG. 1: Schematic diagram of a filament with spatially vary- 
ing spontaneous curvature showing the notation used. 



curvature. In Sec. V we calculate the force-extension 
relationship and identify the various force scales in the 
problem, while in Sec. VI we investigate the response 
of the shape of transverse fluctuations on the stretching 
force. We conclude and discuss a possible extension of 
this work in Sec. VII. 



II. MODEL 

We treat a semiflexible polymer with spontaneous (in- 
trinsic) curvature by simply extending the classical WLC 
model. The filament is described as a plane curve, r(s), 
where s € [0, L] is the arc-length along its contour, sub- 
ject to a local inextensibility constraint so that the total 
contour length L is constant. The constraint is mathe- 
matically expressed by the condition: 



dr(s) 



(i) 



The effective free-energy functional ( "Hamiltonian" ) of 
the polymer depends only on its curvature (bending), 



H [v(s)] 



ds 



dt(s) 



Os 



c(s) 



(2) 



where t(s) = dr(s)/ds is the tangent vector and c(s) 
is the spontaneous curvature at the position s, cf. Figjl] 
This bending energy penalizes conformations with a local 
curvature if it deviates away from c(s). 

The interaction with a stretching force f which pulls 
at the polymer end-points is expressed by adding to Hq 
the work function 



W[v(s)} = -f • [r(L) - r(0)] 



dsf-t(s). (3) 



The stretching response of a filament described by Eqs. 
([2]) and ([3]) is in general analytically intractable for two 
reasons. One is the well-known problem associated with 
the local inextensibility constraint, the other is the non- 
analyticity introduced by the definition of curvature as 
the modulus of dt(s)/ds. However, in the present case of 
weakly-bending filaments in two dimensions the problem 
simplifies drastically and becomes linear. The coordi- 
nate system that we use is shown in Fig. [I] To describe 
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a weakly-bending filament, it is convenient to use the 
parametrization r(s) = [s — X(s), y(s)] which implies 
t(s) = [1 — d s X, t±(s)]. The inextensibility condition can 
be written as t\ + (1 — d s X) 2 = 1 which, in the weakly 
bending approximation "C 1, yields d s X 



1+2 



*i + 



2"_L 
+5 



O(^) and therefore = + ^d a t± + 0{t\). 

In this approximation, the principal normal to the fil- 
ament can be considered perpendicular to its backbone 
(x-direction) and dt/ds ss [0, d s t±]. Introducing the 
signed curvature, c(s), we can express the energy of a 
weakly-bending stretched filament as 



n[y(s)} 



ds 



+ -f f ds 

r Jo 



(4) 



c(s) is positive if the filament bends convexly and nega- 
tive if it bends concavely. 



III. GROUND STATE 

In the absence of thermal fluctuations, the response 
of the filament to the stretching force is determined by 
minimizing the functional H[y(s)] for appropriate bound- 
ary conditions. In this paper, we assume hingcd-hinged 
boundary conditions which can be realized by attaching 
beads at the end-points that can exert forces but not 
moments [26]. As we show below, in the case of strong 
stretching, the specific boundary conditions may become 
irrelevant. Formally, they are expressed by the following 
relations: 



y(0) = y(L) = 
y"(0) = 5(0), y"(L) 



c(L) 



(5) 



(The primes denote derivatives with respect to s.) 
The minimizing configuration, y a (s), obeys an Eulcr- 
Lagrange equation which can be obtained by adding an 
infinitesimal variation, y(s) — y a (s) + £f](s), where e is 



infinitesimal and j](s) is an arbitrary function of s sat- 
isfying »7(0) = rj(L) = and r)"(0) = <q"{L) = and 
requiring the derivative of H(e) as a function of e to van- 
ish. The athermal configuration is the solution to the 
equation 



fy':(s) + KyZ"(s)-Kc"(s) = 



(6) 



with the boundary conditions given by Eq. ([5]). 

We first consider the case of constant curvature: the 
value c(s) = c is maintained along the chain. The last 
term in Eq. ([6| vanishes and the athermal configura- 
tion y a (s) is determined by hyperbolic functions with the 
characteristic length scale I / = \J kJJ. Implementing the 
boundary conditions, we obtain 



y a (s) 



c k 



/sinh(^) 



sinh ( — ) + sinh ( 



(-£_+£)- 
It ' 



-sinh(^) 



(7) 



As expected, in the vanishing-force limit / — > the 
above expression yields a parabola, y a {s) = ^cqs(s — L), 
which is identical, to leading order in the weakly bend- 
ing approximation, to a circular arc of curvature Cq. For 
If <C L, and lf<^s<^L — lf,we get a uniform displace- 
ment, y a (s) = — cq/c//) which implies that the actual 
response comes from the segments of arc-length / f at the 
boundaries. 

For a filament with arbitrary variable spontaneous cur- 
vature, one can expand it in Fourier modes: 



c(s) 



nn 



, cos(qs + <p n ), q = — , n € N . (8) 



Since Eq. ^ is linear, its solution for the general case 
will be the superposition of the modes which correspond 
to the individual harmonics. Taking a particular mode, 
c(s) = c„ cos(<?s + 0), the ground-state stretched config- 
uration becomes: 



y a (s) = 



h (.' 1 1 



(/ + (Z 2 K ) sm h(£) 



,,. . , ,1L, , ,1 (L- 2s). 

2cos(0)smh {-j- ) cosh (- J — cos(qs + • 

2 If 2 If 



i sinh(^) 
l f 



(9) 



r 



Since Eq. ^ is linear, a chain with spontaneous curva- 
ture containing many harmonic modes yields a response 
(y a (s)) which is the superposition of the responses for the 
individual modes. The above expression becomes more 



transparent in the case where c(s) = Cosin(gs), 

cqk sin(gs) 



y a (s) = -- 



f + q*K 



(10) 



In the strong stretching regime, If -C L, any mode yields 
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a uniform response for l/<s<i - It, 

KCOs(qs + (j) 



Va(s) 



f + q 2 K 



(11) 



consisting of a simple reduction in the amplitude. We 
see that the length scale I f defines the distance along the 
chain over which the boundary effects persist (boundary 
layer). 



IV. RESPONSE AT FINITE TEMPERATURE 

Thermal fluctuations are added to the athermal 
ground-state response, y(s) = y a {s) + u(s), and con- 
tribute according to a Boltzmann weight associated with 
the energy given in Eq. Q. As we show below, in the 
linear theory of the weakly-bending approximation, the 
two Gaussian contributions (thermal and athermal) de- 
couple. The elastic energy of thermal fluctuations then 
takes the form, 



HWs)] = lj a d S {u"{ S )f+ f -J^ ds{u'{s)Y 
+k ( L d S u"{s) (y»(s)-c(s))+f [ L dsy' a (s)u\ S ) 



const 



(12) 



Using the Euler-Lagrange equation ([6]), and the bound- 
ary conditions one can easily show that the second 
line of the above equation vanishes. Therefore (u(s)) = 
(u'(s)) = (u"(s)) = 0, where (...) denotes average with 
the weight exp(- PH[u(s)]) , j3 = l/k B T. 

We calculate (u 2 (s)} and ((u'(s)) 2 } by decomposing 
u(s) into Fourier harmonics which respect the prescribed 
boundary conditions: u(s) = ^ n a„ sin(n7rs/L), n £ N. 
This gives: 



r, 2^ 



cos 2 ( - 



KT) 2 + /1 



(13) 



and 



({u(s)) 2 ) = 2k B TL 



OO 



sin 2 (^s) 



(14) 



Note that the spontaneous curvature c does not affect 
these thermal fluctuations. The decoupling of thermal 
and athermal curvature undulations and the results that 
we have obtained so far allow us to calculate various ex- 
perimentally relevant observables such as force-extension 
curves and the average width of transverse fluctuations. 



V. FORCE-EXTENSION RELATIONSHIP 

The average projected length of the filament in the 
direction of the stretching force is given by 



(x(L)} = L 



1 



ds(t 2 As)) 



1 



= L~\j ds (y'a( s )f - \ J ds{(u'(s)) 2 ),{15) 

where we have used the weakly bending approximation. 
There are four different force scales in the general prob- 
lem, 



f = — f = — 

Jcr — r2 ' ■> c — —2 



ft 



K 



h 



1 

K 



(16) 



associated with the corresponding length scales: the per- 
sistence length L p — 2K/ksT 1 the radius of sponta- 
neous curvature 1/co, the wavelength of modulation of 
the spontaneous curvature 2ir/q, and the contour length 
L. The force scale /l is related to finite-size effects in the 
sense that for / 3> /l, or equivalently for If L, the 
filament can be viewed as an effective freely jointed chain 
consisting of "pieces" each of contour length I f [331 HI] • 
In principle, the length scales 1/cq and 1/q are indepen- 
dent (the former is the amplitude of the spontaneous cur- 
vature whereas the latter is the wavelength). However, 
the case where Co > 9 corresponds to very steep un- 
dulations and goes beyond the weakly bending approxi- 
mation in the absence of a strong stretching force. The 
steepest undulations compatible with this approximation 
have Co ~ q. On the other hand, we can be within the 
weakly bending approximation with cq > q if the stretch- 
ing force is strong enough. 

Let us first consider the case of constant spontaneous 
curvature c = cq. In the limit where / 3> Jl, the force- 



extension relationship (151 acquires the simple form, 



W)) = l f 



L 



3/2 



1 (U 

2l / 



1/2 



(17) 



The first term in the right-hand side (rhs) of Eq. (17 1 is 



associated with athermal stretching whereas the second 
term is associated with the stretching of thermal fluctu- 
ations. The crossover force scale \{L/cQ 1 ){f^ 2 / fll 2 ) ~ 
kL p Lcq determines the relative weight of the two con- 
tributions. In the strong stretching regime, the second 
term (which does not depend on the spontaneous cur- 
vature) becomes dominant in agreement with a similar 
result obtained recently in Ref. [19]. If L p « L and 
the spontaneous curvature of the unstretched filament is 
small, Cq" 1 > L, the athermal contribution is almost neg- 
ligible. It becomes significant, however, when the spon- 
taneous curvature is large (the weakly bending approxi- 
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FIG. 2: Force-extension curve (solid line) for a filament with 
Lp/L — 1 and cqL = 12. For comparison, we show the force- 
extension curve of the corresponding filament without spon- 
taneous curvature (dashed line). 



FIG. 3: Force-extension curve (solid line) for a filament with 
Lp/L — 0.05 and CoL p = qL p = 50. For comparison, we 
show the force-extension curve of the corresponding filament 
without spontaneous curvature (dashed line). 



mation still applies to the stretched configuration of the 
filament). This is illustrated in Fig. [2] 

We now consider the case of a filament with sponta- 
neous curvature which varies along its contour. If the 
local spontaneous curvature vanishes at the end-points, 
we can expand it in a sine series. In this case, the force- 
extension relationship takes the form, 



L 



= 1 - 



1 K 7T 



2^2 



l k B T 
2 L 



OO 

E 

n=l 



2 L 2 
1 



E 



K 



(18) 



We see that as the stretching force increases, it succes- 
sively irons out the larger length-scale (softer) modes of 
curvature (spontaneous or thermally excited) [25] . If the 
spontaneous curvature has only one mode and / ^> fa, 
we obtain: 



L 



I- 1 



fqfc 



2 [/, + f? 



J cr \ ' 

7 



(19) 



This result holds irrespective of the value of the local 
spontaneous curvature at the end-points. As expected, 
when the spontaneous-curvature undulations have been 



ironed out by a strong force, the force-extension relation- 
ship is determined by the residual thermal fluctuations. 
The same is true when the persistence length is so small 
that thermal fluctuations wash out the spontaneous cur- 
vature. However, when L p ^> c^ 1 sw q~ x , the sponta- 
neous curvature strongly affects the response to a stretch- 
ing force. We point out that for fa -C f cr -C / ^ f c & fq 
the response will be quite soft despite the nonlinear 
stretching of thermal undulations due to f cr <C / as il- 
lustrated in Fig. [3] If the spatially variable spontaneous 
curvature contains many harmonic modes, our analysis 
still applies: Eqs. ^ and (151 yield a result similar to 
that of Eq. ( 18 ) but with the first sum in the rhs replaced 



by a more complicated term. For fa > /, the value of 
the local spontaneous curvature at the end-points only 
slightly affects the response as shown in Fig. [2] 

We can gain insight into the effect of the spatially mod- 
ulated spontaneous curvature on the strong stretching 
response by considering the effective freely jointed chain 
mentioned above. A two-dimensional freely jointed chain 
consisting of N "pieces" of end-to-end distance 6, has a 
force extension relationship, 



(x(L)) 



= 1 - 



k R T 



(20) 



Nb ' 2/6 ' 
where / 3> ksT/b [53]. This equation can also be written 
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VI. TRANSVERSE FLUCTUATIONS 



0.99 



0.98 



0.97 



0.96 



0.95 




Apart from the force-extension curve, the shape of 
transverse fluctuations can be a useful diagnostic tool 
for the elasticity of a semiflexible filament For a 

fluctuating filament with spontaneous curvature c(s) = 
c sin(gs) (q = wr/X), the average width of transverse 
fluctuations along the polymer contour is given by 



ds((y(s)) 2 ) = 



Kfc 



Lk B T 



if) 



3/2 



(f + f q ) 2 

3arth(*/7)/+3v7 



6/l(/) 5 / 2 



(22) 



where the non-dimensional shorthand / = f / ft is used. 
The first (athermal) term is related to the spontaneous 
curvature and the second term comes from the thermal 
fuctuations. For / S> /l, the latter can be approximated 
by \Lk B T/f. 

The profile of transverse fluctuations has the general 
form: 



FIG. 4: Force-extension curve for a filament at T = with 
spontaneous curvature c(s) = ^sin(^ !L s) (solid line) and 
c(s) = 



2 - 6 cos(^s) (dashed line). 



((y(s)) 2 

L 
L 

+ 2* 



r 2r2 



(f/h) 

G(f/f L ,s/L) 



2 F(f/f L}q ,s/L) 



(23) 



as 



L 



_b_ 
2l f 



J cr \ ' 

7 



(21) 



where the chain is viewed as a succession of N "pieces" of 
arc-length 21 f. In the case of a WLC without spontaneous 
curvature, b = 21 f 23J. In the case of a filament with 
spatially modulated spontaneous curvature, however, the 
"pieces" will be wavy with b < 21 f. Plugging the result 
of Eq. (11) into the athermal part of Eq. 
assuming q~ x <C 21 
f q fc/2[f q + f} 2 



if, we obtain b 
thus recovering Eq. ( 19 1 



(|15|, and 

x(2l f )) « 21 f (I - 



We note that in the case of spatially varying sponta- 
neous curvature, the force-extension response for / 3> /l 
comes mainly from ironing out undulations in the bulk. 
Therefore we do not expect clamped boundary condi- 
tions of small compliance (e.g., with the direction of the 
clamped segments close to that of the pulling force), to 
significantly affect the response. On the other hand, as 
shown in Eq. |7]), for a filament with constant sponta- 
neous curvature and large persistence length in the strong 
stretching regime (/ 3> /l), the response comes mainly 
from the boundary layers and the results will depend on 
the specific boundary conditions. 



where F involves or f q depending on whether the 
spontaneous curvature is constant or oscillating. In 
the absence of spontaneous curvature, the profile of 
transverse fluctuations is determined by the function 
G(f/fL,s/L) oc ((u(s)) 2 ) given in Eq. (fill) which is 



bell-shaped and its height decreases rather quickly with 
increasing force. This also applies when the stretching 
force is so strong that it has ironed out the undula- 
tions due to the spontaneous curvature. Interestingly, 
this behavior persists even for stretching forces such that 
fh ^ fcr <^ / <^ fc ~ fq where, as we have shown, 
the force-extension relationship is strongly affected by 
the spontaneous-curvature contribution. The behavior 
becomes qualitatively different in the case of rather stiff 
filaments (L p ^> L) for / <§; f q ~ f c and f/f q ^> L/L p . 
The athermal term becomes dominant and it decreases 
very slowly with increasing force. In addition, the profile 
is essentially flat, characterized by the constant ampli- 
tude of the sinusoidal undulations. This type of behavior 
is expected also in the case where f c 3> f q and / <C f q 
such that / c // 2 > L/(L p f). In order to get f c > f q , we 
need a filament with very steep undulations of the spon- 
taneous curvature (co ^> q). This would take us out of 
validity of the weakly bending approximation. Although 
we cannot make quantitative predictions for this case, 
we expect the main qualitative features of the predicted 
behavior to persist. 
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VII. CONCLUSIONS AND OUTLOOK 

In the present work, we considered a straightfor- 
ward extension of the WLC model which incorporates 
spontaneous curvature to study the response of a two- 
dimensional weakly bending filament to a pulling force 
applied at its ends. The response of a filament with 
intermediate flexibility (L/L p w 1) and constant spon- 
taneous curvature is affected appreciably only for large 
initial curvature, where the effective stretching modu- 
lus of the chain softens. If the spontaneous curvature is 
inhomogeneous, characterized by undulations of a char- 
acteristic wavelength, we find a range of forces where 
these undulations soften the response despite the non- 
linear flattening out of thermal fluctuations. The pro- 
file of transverse fluctuations is mostly dominated by the 
bell-shaped thermal contribution unless the spontaneous 
curvature undulations are very steep. 

Although our quantitative results are based on the two- 
dimensional weakly bending assumptions, the qualitative 
features are expected to hold more generally as they are 
related to the length and force scales of the problem. As 
we have shown in Sec. V, the softening of the stretch- 
ing response for f L < f cr < / < f c w /, can be 
attributed to the fact that the "pieces" of the effective 
freely jointed chain have an end-to-end distance smaller 
than the corresponding contour length because of their 
spontaneous curvature. One might expect a similar hi- 
erarchy of length and force scales to persist in the case 
of a three-dimensional stretched helical filament with the 
pitch length and the radius playing roles analogous to q~ l 
and Cq , respectively, thus yielding a similar softening. 
Checking the validity of this expectation is an interesting 
direction for future work. 
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